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Abstract. We prove that the vanishing of 11 focal values is not suffi- 
cient to ensure that a plane cubic system has a center. 



1. Introduction 



In 1885 Poincare asked when the differential equation 

/ ^ x + p{x,y) ^_ P{x,y) 
y + q{x,y) ' Q{x,y 

with convergent power series p{x, y) and q{x, y) starting with quadratic 
terms, has stable solutions in the neighborhood of the equilibrium solu- 
tion = (0,0). This means that in such a neighborhood the solutions 
of the equivalent plane autonomous system 

x = y + q{x,y) = Q{x,y) 

y = -X- p{x, y) = -P{x, y) 

are closed curves around (0,0). 

Poincare showed that one can iteratively find a formal power series F = 
x'^ + y'^ + h{x, y) + fi{x, y) + ... such that 



with Sj polynomials in the coefficients of P and Q. If all Sj vanish, and F 
is convergent then is a constant of motion, i.e. its gradient field satisfies 
Pdx+Qdy = 0. Since F starts with x'^+y'^ this shows that close to the origin 
all integral curves are closed and the system is stable. Therefore the sfs are 
called the focal values of Pdx + Qdy. Often also the notation r]2j ■= sj is 
used, and the rji are called Liapunov quantities. 

Poincare also showed, that if an analytic constant of motion exists, the focal 
values must vanish. Later Frommer |Fro34] proved that the systems above 
are stable if and only if all focal values vanish even without the assumption 
of convergence of F. (Frommer's proof contains a gap which can be closed 



Unfortunately it is in general impossible to check this condition for a given 
differential equation because there are infinitely many focal values. In the 
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case where P and Q are polynomials of degree at most d, the sj are polyno- 
mials in finitely many unknowns. Hilbert's Basis Theorem then implies that 
the ideal loo = {si, S2, ■ ■ ■) is finitely generated, i.e there exists an integer 
m := m{d) such that 

Si = S2 = ■ ■ ■ = Sm{d) = Sj = Vj. 

This shows that a finite criterion for stability exists, but due to the indirect 
proof of Hilbert's Basis Theorem no value f or m{d ) is obtained. In fact even 



today only m(2) = 3 is known. Zoladek Zol95j and Christopher |Chr05] 
showed that m(3) > 11. Since the number of variables for d = 2 is six and 
m(2) = 6 — 3 it has been conjectured that for d = 3 with 14 variables one 
has m(3) = 14-3 = 11. 

It is the purpose of this note to prove m(3) > 12. 

The most naive approach to this problem is to calculate a Grobner Basis of 
-^11 = (si, . . . , sii) and prove that si2 hi by the usual ideal membership 
test. Unfortunately this is not feasible, since the sj are very complicated. 
They involve 14 variables and are of weighted degree 2j. For example S5 
has already 5348 terms and takes about 1.5 hours on a Powerbook G4 to 
calculate. The polynomials Sj, j > 6 can not at the moment be determined 
by computer algebra systems. 

Zoladek and Christopher therefore deduce their result geometrically. They 
exhibit a component l^ii C X^o = V{Ioo) that has codimension 11 in the 
space of all possible (P, Q) of degree at most three. Finding a component 
of codimension 12 is not an easy task, and indeed we choose a different 
approach. We prove that there exist a codimension 11 family plane au- 
tonomous system of degree 3 with a focus for which nevertheless the first 1 1 
focal values vanish, but 12th one doesn't. For this we look at the system 

X = y + 3x^ + 8xy + 5y^ + 3x^ + 2bx^y + 20xy^ + 18y^ 

y = -i^x + 21 x^ + 9xy + 22y^ + llx^ + 20x^y + 4xy^ + 3y^) 

and prove that for this system sj = mod 29 for j < 11 while S12 7^ 
mod 29. Checking that furthermore the Jacobian matrix of si, . . . , sn has 
full rank modulo 29 for this system, we can apply a theorem of Schreyer 
jSch96] to show the existence of the desired family of foci over C. From this 
we deduce that S12 In = (si, . . . , sn). If fact we even prove the stronger 
result S12 rad In. 

Since for given a given system one can evaluate the sj using Frommers al- 
gorithm |vBC07] without knowing the complete Polynomials, this approach 
is feasible. 

We found the above system by performing a random search. Heuristically 
each Si vanishes mod 29 for about one of every 29 differential equations 
|vBS05j . So we expect to find an example as above after checking 29^^ ~ 
10^^ random examples. By parametrizing si and S2 we can improve this 
to 29^ ~ 10^^ random examples. Indeed we found the example after about 
8 X 10^^ trials. Using an improved version |vBK09] of the program |vBC05] 
this took 1246 CPU-days. Since this search is easily parallelizable we could 
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do this calculation in about one month by distributing the work to several 
computers. 

We would like to thank the Regionales Rechenzentruni fur Niedersachsen 
(RRZN) and the Institut fiir Systems Engineering, Fachgebiet Simulation for 
providing the necessary CPU time. Also we are grateful to Colin Christopher 
who checked our example using REDUCE |Hea04] . 

2. The Proof 

Notation 2.1. If / C . . . , Xn] is an ideal and Xz = V{I) C is the 
variety over spec Z defined by /, then we denote by Xf^ the fiber of Xz over 
¥p for any prime p. Furthermore we donote by Xq the variety defined by / 
over C. 

Theorem 2.2 (Schreyer). Let I = (/i, . . . , /fc) C Z[xi, . . . , Xn] be an ideal 
and Xz = V{I). If x £ Xf^ is a point with codimTxp^.j; = k then there 
exists an irreducible component Yz C Xz with x £ Yz and Yz <^ Xf^. In 
particular Yq ^ 9 

Proof. This is a special case of a theorem of Schreyer |Sch96| . See also 
IvBELOB] for a proof. □ 




Figure 1 . A variety over spec Z 

Example 2.3. Consider Xz = V{3x) C A^. This variety has two compo- 
nents over Z namely Yz = V{x) and Zz = V{3). Since 3 = is true only in 
F3 we have Zz = Z^^. Furthermore Zq = 0. On the other hand x = is 
possible over all ¥p and Ic / 0- See Figure [T} 

Indeed, if we consider the point x = £ ^Wpi P 7^ 3 then we have that the 
derivative (3x)' = 3 7^ and the tangent space Tq^Xf has codimension 1. 
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Therefore the Theorem apphes and the component 1^ containig x = O^p is 
not contained in Xf^. 

Since 3 • 1 = G F3 we can also consider the point x = 1 G Xf^. Here 
we have (3x)' = 3 = and the tangent space Ti^x^,^ has codimension 0. 
Hence the Theorem does not apply, and indeed the component Zz = Zf^ 
containing x = Ipg is completely contained in Xf^. 



Corollary 2.4. If in the situation of Theorem 2.2 we have a further poly- 
nomial g G . . . , Xn] satisfying g{x) 7^ G Fp then g does not vanish on 
Xr. 



Proof. Assume to the contrary that g vanishes on Xc- By Theorem 2.2 



we 



have a component Yz C Xz with G Iz and Ic 7^ O- Since g vanishes on 
Xc and 7^ is also vanishes on Ic and therefore on Yz and Ip^. But 
this contradicts our assumption g{x) 7^ 0. □ 

Theorem 2.5. m(3) > 12. 



Proof. Use our implementation of Frommers algorithm |vBC07j . |vBC05j . 
|vBK09j or REDUCE [Hea04] to check that the example in the introduction 
satisfies the conditions of Corollary 2.4 □ 
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